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Abstract. The entanglement quantum properties of a spin- 1/2 Ising-Heisenberg 
model on a symmetrical diamond chain were analyzed. Due to the separable 
nature of the Ising-type exchange interactions between neighboring Heisenberg 
dimers, calculation of the entanglement can be performed exactly for each 
individual dimer. Pairwise thermal entanglement was studied in terms of the 
isotropic Ising-Heisenberg model, and analytical expressions for the concurrence 
(as a measure of bipartite entanglement) were obtained. The effects of external 
magnetic field H and next-nearest neighbor interaction Jm between nodal Ising 
sites were considered. The ground-state structure and entanglement properties of 
the system were studied in a wide range of the coupling constant values. Various 
regimes with different values of the ground-state entanglement were revealed, 
depending on the relation between competing interaction strengths. Finally, some 
novel effects, such as the two-peak behavior of concurrence versus temperature 
and coexistence of phases with different values of magnetic entanglement were 
observed. 
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1. Introduction 

During the last two decades low-dimensional magnetic materials with competing 
interactions or geometrical frustration have become an intriguing research object. 
Particularly, these materials exhibit a rich variety of unusual ground states and 
thermal properties, as a result of zero and finite temperature phase transitions [THS]. 
As attractive models among these systems, one should mention the ones, having a 
diamond-chain structure. The latter consists of diamond-shaped topological units 
along the chain (Fig. [T]). It has been observed that the compounds A3Cu3(P04)4 
with A=Ca, Sr [B], Bi4Cu3V20i4 |7j and Cu3(Te03)2Br2 [5] can be nicely modeled 
by the Heisenberg diamond chain. Besides, recent experimental results on the natural 
mineral azurite (Cu3(C03)2(OH)2) tSj showed that Cu^"'' ions of this material form 
a spin- 1/2 diamond chain. Furthermore, the discovery of a plateau at 1/3 of the 
saturation value in the low-temperature magnetization curve [9l[T0] has triggered an 
intensive interest in the magnetic properties of azurite |1 H I12 ) . Azurite falls into 
the class of geometrically frustrated magnets. However, the question of the strength 
and the type of exchange interactions for this natural mineral, despite the long- 
standing interest, is still open. The first diamond spin chain was explored under 
a symmetrical condition Ji — J3 [T3] that predicted magnetization plateaus both at 
1/3 and 1/6 of saturation (14] . The frustrated diamond chain with ferromagnetic 
interactions Ji , J3 < and antiferromagnetic interaction J2 > was also investigated 
theoretically [TS]. Other exchange interactions, like an additional cyclic four-spin [TB] 
and Jm interaction between monomeric units (the so called generalized diamond 
chain) jl7j were considered. Additionally, the importance of an anisotropic exchange 
and Dzyaloshinskii-Moriya interaction [18] or interchain coupling [19] was discussed. 
To sum up, note that the theory predictions for certain values of exchange coupling 
constants within a relatively broad range can fit the experimental results. The 
controversy on these values seems to be cleared up only recently (the latest comparison 
of experimental and theoretical results can be found in Ref. [20]). 

Motivated by the controversies discussed above and the fact that different 
compounds can be described by means of a diamond chain, we shall explore 
systematically the generalized symmetrical spin- 1/2 diamond chain with various 
competing interactions in a magnetic field. Unfortunately, the rigorous theoretical 
treatment of geometrically frustrated quantum Heisenberg models is difficult to fulfil. 
The problem arises due to a non-commutability of spin operators involved in the 
Heisenberg Hamiltonian. This is also a primary cause of a presence of quantum 
fluctuations. Owing to this fact, we will use the recently proposed geometrically 
frustrated Ising-Heisenberg diamond chain model STj-HS]- The latter suggests to 
overcome the mathematical difficulties by introducing the Ising spins at the nodal 
sites and the Heisenberg dimers on the interstitial decorating sites of the diamond 
chain (Fig. [l}. For understanding of the properties of underlying purely quantum 
models it is required to obtain an analytic expression for all thermodynamic functions 
of the model. Note that some exactly solvable models with Ising and Heisenberg bonds 
can also provide satisfactory quantitative picture |24| . 

In the present paper we shall mainly deal with the quantum entanglement 
properties of the spin-1/2 Ising-Heisenberg model on a generalized symmetrical 
diamond chain. It is well-known, that the entanglement is a generic feature of 
quantum correlations in systems, that cannot be quantified classically j251l26j . It 
provides a new perspective for understanding the quantum phase transitions (QPTs) 
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and collective phenomena in many-body and condensed matter physics. This problem, 
which has been under scrutiny for nearly two decades, has recently attracted much 
attention pTMTj . A new line research points to a relation between the entanglement 
of a many-particle system and the existence of the QPTs and scaling 32,33J. On 
the other hand, the study of entanglement in solid state physics is of a great 
relevance to the area of Quantum Information and Quantum Computation, since 
many proposals of quantum chips are solid state-based. Although it was believed 
that the entanglement should not manifest itself in macroscopic objects (because of 
a large number of constituents interacting with the surroundings that induce the 
decoherence phenomena), it was theoretically demonstrated that entangled states can 
exist in solids at finite temperatures. This kind of entanglement is referred to in 
literature as "the thermal entanglement" ^7\. And afterwards a few experimental 
evidences have been reported for low-dimensional spin systems |34| . confirming the 
presence of entanglement in solid state materials. 

Returning to the spin- 1/2 Ising-Heisenberg model on a diamond chain we remark 
that the nodal Ising spins represent a barrier for quantum fluctuations. On the other 
hand, taking into account that each Heisenberg dimer interacts with its neighboring 
dimer through the Ising-type, i.e. classical exchange interaction, we find that the states 
of two adjacent dimers become separable (disentangled) Thus, we can calculate 
the concurrence (the measure of pairwise entanglement |35)), which characterizes 
quantum features of the system, for each dimer separately. The main objective of 
the paper is to reveal different regimes of the symmetrical diamond chain and to 
analyze new quantum effects (such as double peak behavior in the concurrence versus 
temperature curves, existence of magnetic entanglement |27j of two different values). 

The rest of the paper is organized as follows: we start in Sec. [5] by obtaining 
concurrence as a measure of entanglement of the spin- 1/2 Ising-Heisenberg model 
on a generalized symmetrical diamond chain. The ground-state structure and the 
entanglement features of the ideal diamond chain ( Ji = J3, J,„ — 0) are discussed in 
Sec. 12.11 The following section contains similar results with the incorporation of 
interaction. Some comments and concluding remarks are given in Sec. HI 

2. Concurrence and thermal entanglement of the spin-1/2 
Ising-Heisenberg model on a generalized symmetrical diamond 

We consider the spin-i Ising-Heisenberg model on a generalized symmetrical diamond 
chain (Ji = J3 = J), which consists of monomeric and dimeric sites (empty and 
full circles in Fig. [U respectively). Within the proposed Ising-Heisenberg model, 
the monomeric (nodal) sites are occupied by Ising spins, while the dimeric sites by 
Heisenberg-type spins. The Hamiltonian can be written as follows: 
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where the summations run over clusters (Fig. [T|), T-Lk represents the Hamiltonian of 
the k — th cluster, Sfc = (S^^, Sf, S^) denotes the Heisenberg spin-i operator, fik 
is the Ising spin. Considering, that each Ising spin belongs simultaneously to two 
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Figure 1. A cross-section of a generalized symmetrical diamond chain 
(k labels the number of the cluster). The empty (monomeric units) and 
full circles (dimeric units) denote lattice positions of the Heisenberg and 
Ising spins (within the proposed Ising-Heisenberg model), respectively. 
Solid lines schematically reproduce the Heisenberg J2 interactions between 
dimeric units, while the broken ones label the Ising-type (nearest-neighbor 
Ji, J3 and next-nearest neighbor J,„) interactions. 



clusters, we have taken a 1/2 factor for the Ising spins in the last term of (1), which 
incorporates the effects of external magnetic field. J, J2, Jm > corresponds to the 
antiferromagnetic couplings. The system will be strongly frustrated due to the chain's 
geometry and existence of competing interactions J, J2 and Jm- When J,„ = we deal 
with the so called ideal diamond chain [21 . Before introducing the calculations and 
discussion we would like to emphasize the fact which was already discussed in Sec. [T] 
the states of two neighboring Heisenberg dimers (with interaction J2) are separable 
(disentangled), because of a classical character of the coupling between them (by means 
of the Ising spin). Hence we can calculate the entanglement for each of the dimers 
individually. Note, that a different approach of a mean-field-like treatment, based on 
the Gibbs-Bogoliubov inequality was used in Ref. [36_ , where all the couplings between 
the diamond chain sites were chosen to be of a quantum (Heisenberg-type) character. 

Here we use concurrence (as a measure of bipartite entanglement |35j ) of the 
Heisenberg dimers by tracing out the Ising spins in each cluster. For the construction 
of eigenvectors of each cluster we will take into account that Hk possesses a symmetry 
corresponding to the permutations /ifc -r-?> /ife+i and {/ifc O ^k+i',Sk^ O Sfca}- 
Besides, the Hilbert space of the cluster Hciuster can be presented as Hciuster = 
T-ik (8) Hdimer ® "Hfc+i, whcrc Hki, Hdimer, T~Lk2 dcnotcs the Hilbert spaces of ^fc, 
Heisneberg dimer and /ifc+i respectively. We obtain the following eigenvectors, due to 
the symmetries and Hilbert space structure (hereafter, the letter k labels the number 
of the cluster): 

= ^(1 tfci ® ik/j + Ufei «) tfe,)) ® I tfctfe+i); 
1^2) = ^(1 tfci ® ifcj + Ufei ® tfe2)) ® (I tfelfe+i) + I Ifctfc+i)); 
ItZ-s) = tfci ® ifc2) + Ufci «) t/c^)) «) (I tfeifc+i) - I ifctfc+i)); 
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IV-fcia) = I ifei «)ife2) «) I tfctfe+l); 
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I tfci tfca) «) I tfetfe+l); 



(2) 



and the corresponding eigenvalues: 
1 



Ei = - {-2H + Jm + J2) ; E2=E3 = 
E4 - \{Jrn + J2 + 2H); E5 - h-2H + J,„ - 3J2); 



Ee — Ej 

Eg ^ 



4^ 

— — ; -; Es — -{2H + J,n ~ 3J2); 



3H J„ 



J', E] n — E^ 1 — —H — 



(3) 



1 1 

£^12 = - Y + 4 + J2 - 4J) ; Ei3 = — + - (J„ + J2 - 4J) ; 
Eli = £15 = J,» + J2 + Ei6 = ^+ I + J. 



We study concurrence C{p), to quantify pairwise entanglement |35) . defined as 

C(p) = maa:{Ai - A2 - A3 - A4, 0}, (4) 

where A^'s are the square roots of the eigenvalues of the corresponding operator for 
the density matrix 

/5 = P12 (cr^ «) erf )p:^2(c^l erf) (5) 

in descending order. Since we consider pairwise entanglement, we should use the 
reduced density matrix pi2, by tracing out two (of four) spins of the cluster. The 
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reduced density matrix pi2 is defined as [37] 

Pi2 = ^{a\p\a). (6) 

a 

In this equation \a) denotes basis vectors of the Hilbert space associated with the 
system, with respect to which the density matrix is reduced. The summation runs 
over all these basis vectors. Since in our case we make reduction with respect to two 

spins, |a) ={|U),Ut),in),ltt)}. 

It is obvious that the only entangled pair is formed by the Heisenberg spins. 
Other pairs are disentangled (separable) because of the classical (diagonal) character 
of the Ising-type interaction between them. Hence we will be interested in the reduced 
density matrix, constructed by tracing out two Ising-type spins pL^ and pk+i- In other 
words, pfcj2 = Tr{^j.^^j.^j}pfc and the full density matrix pk of the k — th cluster is 
defined as (here and further Boltzmann constant is set to be fc^ = 1) 

1 

i—l 

where is the partition function: 

6H + .7m+4J+.72 / -fJ + Jm+S J 3H+.J^+2J 

Zk = Trpk = e ^ ( 2e 5T + 2e ^ + 

2e 2T + 2e 2t _|_ g t -|_ e t -f (g) 

2H + J Jf + 2.7 ff+.7 + .72 2H + .7+J2 3H \ 

e T + e T +e t -|_e t + e t + l\ , 

Using the definition ([6]), the basis vectors \a) = { \ H), \ |t), | t|), | tt)} we construct 
the reduced density matrix pk-^^ of the k — th cluster: 
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where 



iH + Jm-J2 -277+ Jm -4.7+J2 -677 + +4J+ J2 

u = 2e + e + e -it 

6H + .7m+4.7+ J2 / H + .7Tn+2J H + 2.7 

V — e 4T ( 2e 2T + e T -\- I 



w 



1 / J2 \ 2H + .7„,+J2 / 77+ \ 

-[e—+lje IT i^2e^T^ +e"'^ + lj , (10) 

1 / 72 \ 277 + .7„ + J2 / 77+ J„. „ /rp \ 

y= - - (^e~ - 1 j e ^ (^2e^^ + e"'^ + 1 j . 

The density matrix pki2 in Eq. (0) has a form of a so called X-state [35], since 
the Hamiltonian 'Hfc is translationary invariant with a symmetry [S'^i'Hfe] — 
(5^ = l/2{pl + pI_^_^) + Sl^ + SI/) [^. The concurrence C{p) of such an X-state 
density matrix has the following form |40| : 

C(p) = |max(|y|- V^,0). (11) 

We note here that the reduced density matrix p^^^ of any pair of spins, different 
from the Heisenberg dimer has no non-diagonal elements, responsible for the quantum 
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Figure 2. Concurrence C(p) versus temperature T for J2 = 1, Jm = 0, 
H = 0, and different values of J. 



correlations, i.e. entanglement (see Thus we conclude, that there is no 

entanglement between a pair of spins which contains at least one Ising spin. 

In Eq. ([3|), one finds a set of states with maximum value of entanglement, for 
which the Heisenberg dimer is in a singlet or a triplet state {ip/s with i = 1, 8). As 
for the rest of the states {ip/s with i — 9, 16) the Heisenberg dimer is in a separable 
state and therefore these ip/s are non-entangled ones. 

2.1. Ideal diamond chain 

In this section we proceed to the investigation of entanglement features of a dimeric 
unit of an ideal diamond chain {Jm = 0). First, we study the behavior of C{p) 
aX H — 0. We will discuss here three regimes, depending on the value of J — J2: 
J — J2 > 0, J— J2 < and J — J2 = 0. In the first case, as one finds out from (3), the 
ground state contains two-fold degenerate states ^^12 and tpi^. Since these states are 
factorable, the corresponding dependency curve of C(p) from temperature T starts 
from C{p) = (Fig. [2]). Furthermore, the entanglement can be invoked by increasing 
the temperature (for values of J — J2 close to 0). This happens since the contribution 
of entangled states in the mixture pk increases with the growth of temperature T. 
The local maximum, appearing here arises due to the optimal thermal mixing of all 
eigenstates in the system. This maximum becomes narrower and smaller and gradually 
vanishes by increasing J — J2 . But the value of J — J2 corresponding to disappearing 
of C{p) also depends on the value of J2 (e.g. for J2 = 1, J — J2 « 0.2). The latter 
becomes obvious, if one takes into account that J2, being the coupling constant of 
the Heisenberg type interaction between dimeric units, is responsible for the strength 
of quantum correlations between Heisenberg spins. We would like to emphasize here 
that in the case J — J2 > the system exhibits weak (0 < J2 < J) or no frustration 
(J2 < 0). 

In the second case, when J—J2 < 0, the system will obviously manifest more of its 
quantum nature. Firstly, the dependency curve of C{p) from temperature starts from 
C{p) ^ I at T = (Fig. [2]). This happens due to the fact that at zero temperature the 
maximum entangled states V'S: V'6j V'7 and tp& form four- fold degenerate ground state 
with the value of C(p) = 1 for the corresponding reduced density matrix pki2- When 
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Figure 3. Sudden-death temperature Td corresponding to the vanishing or 
arising of entanglement at zero magnetic H versus ration parameter J2/ J 
(J = 0.5). 



the temperature is increased, the concurrence gradually disappears because of the 
thermal mixing with other states of the system (including the factorable ones). The 
sudden-death temperature Td, corresponding to the dying out of quantum correlations 
in the system can be found through the equation C{p) = 0. It has the following form: 

x--^ {x^ + lf ^2\x''' (12) 

where x = e^^"^. The solution can be presented in the form Td = J /log a (when 
J — J2 < 0), where a depends on the ratio parameter J2/J. Increasing this 
ratio, a decreases, but the linear dependence on J remains (e.g. when J2/J — 2, 
a=i(3 + VT7)). 

Finally, the case J — J2 = can be regarded as a boundary case in the following 
sense. Here the ground state is six-fold degenerate, containing additionally ■0i2 and 
tpi3, besides ■05, ipe, ipy and ips (in other words all the states as in previous two cases). 
Since the "012 and ipis are factorable, this leads to lower entanglement of the ground 
state's reduced matrix, that is C{p) — 1/3 (Fig. [5]). Moreover, the above discussed 
sudden-death temperature Td is lower, than in the case J — J2 < (although again 
Td = J I log a with a = 2 Vs). 

On the other hand, as it can be seen from Fig. [21 there are two sudden-death 
temperatures in the case J — J2 > (corresponding to arising and vanishing of 
entanglement) |41j . The dependence of Td on the ratio parameter J2/J is shown 
in Fig. [31 In the area < J2/J < 1, there are two sudden-death temperatures (as 
mentioned above), while for the values J2/ J ^ 1, the dependence is a linear one. 

Our further investigation concerns the effects of the magnetic field H. 
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Figure 4. (a) Sudden-death temperature Td corresponding to the vanishing 
or arising of entanglement versus magnetic field H for J = 2 and different 
values of J2; (b) Concurrence C(p) versus temperature T for J = 2 and 
different values of J2 and magnetic field H. 



Firstly we will discuss how magnetic field affects the above introduced sudden- 
death temperature T^,. While increasing H, Td increases too, but it always remains 
lower than J2/log3 (in other words liuiH^ooTd = J2/log3). Another interesting fact 
is that magnetic field gives a rise to more than two sudden-death temperatures in the 
case J—J2 > [Fig.|4](a)] and on the dependence of C{p) from temperature T one finds 
two peaks separated by an area of a zero entanglement [Fig.|4](b)]. With increasing H 
the smaller of aforementioned peaks tears apart from C(p) = 0, starts merging with 
the bigger one and eventually disappears. An effect of this kind has not been reported 
yet, to the best of our knowledge. Although a similar two-peak behavior of concurrence 
was found in dissipative the Lipkin-Mcshkov-Glick model versus magnetic field |42| . 
However, when T — >■ 0, C{p) remains finite and becomes zero only at absolute zero 
temperature T — [i.e. there can be not more than three sudden-death temperatures 
corresponding to disappearing or arising of thermal entanglement, as it can be also 
seen form Fig. |4] (a)]. In other words in the area of low temperatures the behavior of 
concurrence is smooth, in contrast with the case when magnetic field is absent. 

Now, we concentrate on the dependence of C{p) on magnetic field. Because 
of the above introduced ground state structure, the dependency curve of C{p) from 
magnetic field at zero temperature has a dip at H ^ with C(p) = 1/3 for J— J2 = 0. 
There is no dip if J— J2 < (Fig. [5]). When Ising-type interaction is stronger than the 
Heisenberg one ( J— J2 > 0), one does not find a magnetic entanglement. Furthermore, 
magnetic entanglement is of a higher value than that at zero magnetic field in the case 
J— J2 = 0. This happens due to the fact that ground state here is two-fold degenerated 
and contains ip^ and ipi2 with the value C{p) — 1/2 for the corresponding reduced 
density matrix. C{p) becomes zero for the case J — J2 < at the values of H, 
corresponding to saturation field, that is when the non-entangled state tttt (in the 
area H > 0) or IIH (in the area H < 0) becomes the ground state. One can find the 
described values of H from the conditions Eg = £'5 and E'lg ~ Eg, giving i7+ = J + J2 
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Figure 5. Concurrence C(p) versus magnetic field H for different values 
of temperature, J2 and J. 
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Figure 6. Concurrence C(p) versus magnetic field H and temperature T 
for (a) J2 = 2 and J = 2; (b) J2 = 1.7 and Ji = J3 = J = 2. 



and H~ = — J — J2, respectively. Thermal effects smoothes the step-like behavior of 
concurrence in the case when J — J2 > and induces thermal entanglement when 
J — J2 > (see Fig. [5]). The further increase of temperature causes the quantum 
correlations eventually dying out for the both cases. 

Summarizing, in Fig. IH] we also plot three-dimensional dependencies of the 
concurrence C(p) versus temperature T and magnetic field H. 



3. Incorporation of J™ interaction 

In this section we will study the effects of the next-nearest neighbor interaction 
between the Ising spins of the cluster, using the full expression for (10) and pT|) . We 
will start with the discussion of the ground state structure for the case H — and 
J — J2 > 0. It turns out that here one can distinguish two regimes. First, when 
< Jm < 2(J — J2), frustrated ground state contains two-fold degenerate -012 and 
tpi3 and thus the dependency curve of C{p) from temperature starts up at C{p) = 
. However, the thermal effects can cause the thermal entanglement for the values 
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Figure 7. Sudden-death temperature Td corresponding to the vanishing 
or arising of entanglement versus Jm for H = and different values of J2 
and J — 2. 



of Jm, close to 2(J — J2) (but remaining Jm < 2(J — J2)). One finds two sudden- 
death temperatures on the dependency curve of C{p) from T (Fig. [7]). This effect can 
be understood from the following discussion. The ground state consists of four-fold 
degenerate states ^e, -07, "012 and ^13, for the case Jm = 2( J — J2). Although this 
mixture contains maximum entangled states "06 and ipj^ the corresponding density 
matrix for this ground state gives C{p) = 0. By increasing temperature, one obtains 
the thermal mixing of states which leads to a higher contribution of entangled states. 
This contribution, however, becomes less, when the values of J,„ are considerably 
higher than 2( J — J2). Thus, with increasing the difference of Jm and 2{J — J2), the 
local maximum becomes narrower and eventually disappears. 

In the opposite case, when J™ > 2( J— J2), the frustrated ground state is two-fold 
degenerate, but with ■06 and 07e, hence above mentioned curve of C{p) starts from 
C(p) = 1. We find only one sudden-death temperature here, which increases with the 
growth of Jm (Fig.[7|). In other words, the qualitative picture remains the same as for 
the case Jm = 0. 

C{p) is of a maximum value {C{p) = 1) at zero magnetic field and zero 
temperature, regardless of the Jm for a dominant Heisenberg interaction ( J — J2 < 0). 

Concluding the discussion of zero magnetic field properties in the case Jm ^ 0, we 
note, that when J,,,, < (ferromagnetic coupling), the absolute value of J,„ does not 
interfere with the ground state properties of the system (it will be two-fold degenerate 
ipi2 and V13, ii J > J2 or V's and -08, if J < J2)- 

Here then, we will discuss the regime J—J2 > introducing effects of the magnetic 
field H. We differentiate two subcases. First one, when Jm < 2(J — J2), one does 
not find magnetic entanglement in the system, since increasing the absolute value of 
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Figure 8. Concurrence C(p) versus magnetic field H and temperature T 
for J2 = 1 and J = 1.5 and Jm = 1-5. 



magnetic field H, we obtain a sequence of separable states (e.g. V12 — ^ (V'lo + ^11) — ^ 
^/jg or ?/)i2 — >■ V'9 for H > Q). Here and further by {psii + we will mean two-fold 
degenerate states. 

Meanwhile, when Jm > 2{J — J2), the aforementioned sequence of states starts 
from {^Q + "07) with maximum value of C{p) = 1, therefore we obtain magnetic 
entanglement (Fig. |8|). One can introduce here critical values of magnetic field iJ+ 
and H~ , corresponding to vanishing of magnetic entanglement. In contrary with the 
case Jm = 0, does not coincide with the saturation fields (see Sec. 12. ip . 

On the one hand, we have the ground state transitions {ipe + ip'j) — > ipi2 — > 
(V'lo + V'li) ~^ '09 (for -ff > 0) if the value of magnetic field, corresponding to the 
intersection of energies and £^12 is lower than that of Eiq and £^12 [Fig.|9](a)]. This 
condition gives: Jm < 2J — J2. Thus the corresponding critical values of magnetic 
field can be found from Eq = E12 with — 2 J2 — 2 J + (obviously, H~ ~ —H^, 
from equation Ee = E13). On the other hand, when J„i > 2 J— J2, we have the ground 
state transitions (-06 +^7) — ^ (V'lo + 0ii) — ^ V'Q [Fig.[9](b)]. Corresponding i/+ ~ J2, 
found from Eq = Eiq [H^ — — J2, from E^, — En). The ground state transition 
(•06 + "07) ~^ V'g can not occur, since the corresponding condition is inconsistent with 
Jm > 2 J — J2. 

Following the same technique as in previous paragraph (we will not stop on 
detailed phase structure), we obtain the following regimes for J — J2 < 0: = 
H+ = J + J2 {H^ = in the case Jm < J2 - J and H+ = 2J2 - Jm if 

J2 > Jm > J2 ^ J , and finally, when Jm > J2 one finds — J2. 

The special (boundary) case J — J2 = is also of interest, since one can observe 
here magnetic entanglement of different values (C(p) — 1 and C{p) = 1/2) fFig. [T0|) . 
whereas in the case Jm — 0, these two regimes cannot coexist for a fixed values of J 
and J2. This situation arises only for < Jm < J, when one finds the sequence of 
states {ipe + 0?) — > (05 +"012) — ^(factorable state) (for H > 0). In other words, at the 
values of magnetic field H — ±Jm (found from conditions Eq — E12 for iJ > and 
Ee = Es for H < 0) the states with different values of magnetic entanglement coexist. 

As for the sudden-death temperature corresponding to the disappearing or arising 
of entanglement at non-zero magnetic field, one finds a similar behavior as in the case 
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Figure 9. Eigenvalues Ea, Eio and E12 versus magnetic field H for (a) 
J™ < 2J - J2; (b) J^<2J- J2. 




Figure 10. Concurrence C{p) versus magnetic field H and temperature T 
for J2 — 2 and J = 2 and Jm — 1.5. 



Jm = 0, i.e. here again we find up to three sudden-death temperatures (as in Fig.|4]), 
with a two peak behavior on the dependency of C(p) on temperature. 

4. Conclusion 

In this paper we have studied the thermal entanglement of a spin-1/2 Ising-Heisenberg 
model on a symmetrical diamond chain, which has been proposed to understand a 
frustrated magnetism of the series of compounds, like A3Cu3(P04)4 with A=Ca, Sr, 
Bi4Cu3V20i4, Cu3(Te03)2Br2 and Cu3(C03)2(OH)2. We have studied the phase 
structure and entanglement properties of the system in a wide range of Ising-type 
interaction constants Ji = J3 = J, Jm and Heisenberg-type J2, considering that 
diamond chain structure describes a broad class of materials (within different values 
of exchange interaction parameters) and that the exact value of coupling constants for 
azurite (Cu3(C03)2(OH)2) is still under scrutinizing question. Taking into account the 
classical and hence separable character of Ising-type interactions which are coupling 
adjacent Heisenberg dimers, we have calculated the entanglement of each of these 



Thermal Entanglement of a Spin-1/2 Ising-Heisenberg Model on a Symmetrical Diamond ChainXA 

dimers separately. We have used the concurrence for quantifying the amount of 
entanglement between two Heisenberg-type spins, by tracing out Ising-type ones from 
the density matrix of the diamond-shaped cluster (the only entangled pair here is 
the Heisenberg dimer). The incorporation of next-nearest neighbor interaction Jm 
has also been investigated (generalized diamond chain) and the effects of external 
magnetic field have been invoked. 

We have revealed a number of regimes with distinct ground state structure 
and qualitatively different thermodynamic behavior, depending on the relations 
between J, J2 and Jm and values of magnetic field H . We found that in general 
for a dominant Heisenberg-type interaction (J2 > J) the system's ground state is 
maximally entangled, but increasing the temperature, pure quantum correlations 
eventually disappears. On the other hand, for a dominant Ising-type interaction 
(J > J2) the ground-state is non-entangled, whether the temperature gives rise to 
thermal entanglement. In the latter case one does not find magnetic entanglement 
at the absolute zero temperature (the system behaves as a classical one). However, 
magnetic field can lead to another, yet not described effect of two-peak behavior 
of concurrence C[p) versus temperature with three sudden-death temperatures (one 
of them corresponding to reappearance of concurrence and the other two to its 
disappearing). These two peaks are separated by an area of a zero entanglement, 
which becomes narrower, with the growth of the magnetic field and aforementioned 
peaks merge into each other. Another novel effect was indicated for a boundary case 
J = J2 when < Jm < J. Specifically, two states with different values of magnetic 
entanglement coexist for the value of magnetic field H — ± J^. One finds a step like 
behavior of concurrence versus magnetic field H, with plateaus at the value 1/2 and 1. 
In other words, the presence of competing interactions in the system and geometrical 
structure of the chain, each leading to a frustration, makes the phase structure of the 
system richer and gives rise to an interesting physical behavior. Finally, the adopted 
model guaranties an experimental realization for suitable theoretical treatment. Our 
results will be useful for further experimental detection of entanglement in the diamond 
chain structured macroscopic samples by means of entanglement witnesses (e.g. built 
from measurements of magnetic susceptibility [13]). 
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